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ABSTRACT 

Species  densities  of  Nitrogen  plasma  are  investigated  using  MATLAB  simulation  technique.  In  this 
simulation,  a  total  of  nine  reactions  with  known  reaction  rate  coefficients  are  used.  From  these  reactions  we  get 
continuity  equations  which  are  simulated  based  on  the  steady  state  rate  balance  equations  (non-linear  equations ) 
derived  from  the  continuity  equations.  The  whole  simulation  is  performed  on  MATLAB  application  software  using  f 
solve  function.  Then  the  species  densities  at  temperature  300  K  by  the  f  solve  values  are  determined.  Electron  density 
has  been  estimated  from  this  simulation  and  is  in  the  order  10 18  cm'3.  The  densities  of  different  Reactive  Nitrogen 
Species  are  determined  from  pressure  1  torr  to  760  torr  are  tabulated  and  discussed. 
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1.  INTRODUCTION 

Plasma  is  a  quasi  neutral  ionized  gas  roughly  equal  number  of  positive  (ions)  and  negative  (electrons)  ions. 
Like  gas,  it  have  no  fixed  shape  or  volume  and  less  dense  then  solid  or  liquid.  Unlike  gas,  it  has  equal  number 
of  ions  of  opposite  charge.  Plasma  can  be  found  in  two  ways,  naturally  produced  plasma  in  sun,  in  various 
stars  or  in  our  earth’s  atmosphere  and  plasma  created  in  the  laboratory.  Various  methods  can  be  used  to  produce 
plasma  in  laboratory  but  one  of  the  most  effective  and  efficient  way  is  by  dielectric  barrier  discharge  (DBD). 
Di  electric  barrier  discharge  (DBD)  is  a  kind  of  gas  discharge  developed  in  1857  by  Werner  von  Siemens  as  a  process 
for  Ozone  production.  For  its  irrelatively  low  cost,  low  operation  temperature  and  of  high  pressure  it  has 
numerous  numbers  of  applications  like  de-pollution  of  a  confined  place,  biological  applications  and  material 
surface  processing  applications.  It  is  proved  that  DBD  can  be  an  effective  source  of  different  reactive  species 
such  as  03,  NO,  N,  O,  OH  etc.  N  (4S),  N+  and  free  electrons  are  the  main  species  with  energy  conservation 
non-equilibrium  excitation  in  the  post-discharge  period  [1],  In  the  inelastic  collision  there  energy  and  momentum 
is  not  the  same  as  it  was  before  collisions.  Energy  transferred  from  one  species  to  another  during  the  collisions 
between  various  species  creates  excitations.  This  collision  sometime  generates  new  species.  The  generation  of 
new  species  mainly  depends  on  the  kinetics  of  the  various  species.  In  the  real  world  it’s  hard  to  generate  a  specific 
species  of  plasma.  Before  generation  of  specific  plasma,  simulation  is  to  be  performing  under  various 
conditions.  This  simulation  is  time  and  cost  efficient.  In  addition  it  also  gives  us  the  idea  of  instrument  it  required. 
Simulation  can  be  modelled  on  various  ways.  Kinetic  modelling  helps  to  simulate  a  system  that  contains 
various  elements  (such  as  gas,  radicals,  ions,  electron  etc.).  Kinetic  model  is  based  on  the  kinetics  of  the 
elements  in  plasma.  In  this  paper  a  simulation  based  on  some  reaction  in  Nitrogen  plasma  is  performed  in  a  close 
system  at  300  K  to  get  different  species  densities. 
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2.  PLASMA  MODELLING 


During  different  applications  of  plasma,  the  plasma  properties  can  be  divided  into  the  following  subgroups  like 
physical  processes,  gas  phase  kinetics  and  surface  interactions  [2],  For  simulations  a  particular  plasma  process  needs  to 
deal  with  the  above  processes.  Depending  upon  the  methods  of  simulation  on  theoretical  approach  we  see  different  types  of 
simulation  techniques  like  MATHLAB  f  solve,  Monti-Carlo  Simulation  etc.  For  determining  the  plasma  properties 
Numerical  Modelling  is  computer  based  method.  For  structural  and  fluid  investigation  this  method  is  useful  by  the 
application  of  laws  of  mechanics  depending  on  the  number  of  pre -prepared  code  which  can  be  used  for  linear  of  non-linear 
analysis.  In  order  to  reproduce  the  plasma  properties  using  the  effect  of  mechanics  and  to  predict  the  reaction  of  various 
plasma  variables  numerical  modelling  is  a  useful  tool.  Numerical  models  are  also  a  useful  tool  for  modelling  plasma 
systems  that  are  too  complex  for  analytical  modelling  methods  to  be  used  [3,  4], 

2.1  Kinetic  Modelling 

To  describe  plasma,  the  kinetic  model  is  a  fundamental  way,  which  produce  a  distribution  function/ (x,  v,  t), 
where  the  independent  variables  fand  v  are  position  and  velocity,  respectively  [5],  By  solving  the  Boltzmann  equation,  a 
kinetic  description  is  achieved.  When  the  description  of  long-range  Coulomb  interaction  is  needed,  by  the  V  solve  equation 
which  contains  self-consistent  collective  electromagnetic  field,  or  by  Fokker-Plank  equation,  where  approximation  have 
been  used  to  drive  manageable  collision  terms.  Using  Maxwell’s  equations  the  charges  and  currents  produced  by  the  self- 
consistent  distribution  functions  can  be  determined. 


We  get  the  species  velocity  distribution  functions  from  the  solution  variables  in  the  kinetic  models.  It  is  assumed 
for  fluid  models  that  the  distribution  function  is  in  equilibrium.  Then  the  Boltzmann  equation  can  be  replaced  by  its 
corresponding  conservation  of  mass,  momentum  and  energy  [6], 


The  plasma  modelling  in  a  continuum  regime  starts  from  the  Boltzmann  equations  [6] 


Sfi 

dt 


+  ct.  V/  +  F.  Vc 


coll 


(1) 


Where  the  index  i,f,c,F,  V  and  Vc indicate  single  species,  distribution  function,  the  charged  particle  velocity,  the 
force,  the  gradient  operator  in  physical  space  and  the  gradient  operator  in  velocity  space. 


The  subscript  ‘coll/  refers  the  change  of  velocity  distribution  due  to  collision.  F  represents  the  Lorentz  force  on 
the  charged  particles.  A  detailed  discussion  on  the  theory  of  Boltzmann  equation  and  its  relation  for  the  fluid  models  is 
presented  in  [6,  7],  The  index  ‘i’  indicate  for  every  type  of  species.  To  get  self-consistent  solution  of  the  above  we  need  to 
solve  for  the  electric  field  E  by  means  of  Poisson’s  equation 

V.  E  =  —  (2) 

£0 


2.2  Reaction  and  Rate  Balance  Equations 

The  considered  Nitrogen  plasma  reactions  and  there  rate  coefficients  at  300  K  are  tabulated  in  table  1 . 


Table  1:  Reaction  &  Rate  Coefficient  at  300  K 


No. 

Reaction 

Reaction  Rate,  k 

Reference 

1 

N?(A)  +  tf2(aO  -»  A4+  +  e 

4.0  x  10~12 

[8,1] 

2 

JV7(aO  +  JV7(aO  ->  JV4+  +  e 

1.0  x  IQ"11 

[8,1] 
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Table  1:  Contd., 

3 

N?+  +  N2(A)  -»  JV3+  +  N 

3.0  X  10~10 

[9,1] 

4 

Nf  +  N  ->  N?+  +  N2 

6.6  X  10-11 

[8] 

5 

N2(A)  +  N  ->  N2  +  N 

2  X  10“12 

[8,1] 

6 

N2(A)  +  N  ->  n2  +n(  2p ) 

4  x  10-11 

[8,1] 

7 

N?{A)  +  JV2  ->  2N? 

3.0  x  10~16 

[8,1] 

8 

Al2+  +  N  ->  N+  +  JV2 

7.2  x  10~13 

[10] 

9 

N(  2P)  +  N2  ->  N2  +  N 

6  x  10“14 

[10,13] 

All  the  reactions  are  considered  in  the  closed  system  having  at  temperature  300  K.  The  reactions  are  assembled 
into  a  model  for  intrinsic  chemical  kinetics  for  the  concentration  of  these  ten  species  [11].  If  the  reactor  is  assumed  kinetics 
dominated  (well  mixed)-when  there  action  time  is  much  shorter  than  the  am  bipolar  diffusion  time,  then  each  of  the 
elementary  reactions  in  Table- 1  contributes  mass  action  law  terms  to  the  kinetics  model  with  rate  constants,  with  number 
densities  defined  as:  Al^is  the  density  of  N2 ;  NNz(A)the  density  of  1V2(.i4);  Al^^-jthe  density  of  N2(a ');  NNthe  density  of 

N;  AfW(.  2pjthe  density  of  N(  2P );  Afw+the  density  of  N+ ;  NN+the  density  of  N2;  NN+is  the  density  of  N3;  Afw+the  density 
of  Nf  and  Ne  the  density  of  electron  [11]. 

The  species  density  is  a  function  of  time.  The  continuity  equation  is  needed  to  simulate  the  time  dependent 
behaviour  of  the  Nitrogen  plasma  species  densities.  Now  we  can  write  the  continuity  equations  from  the  reactions  and  rate 
coefficients  that  shown  in  table  1  as 


dNN2 

dt 


—  k4  Nn+NNz  +  k5NNlNN  +  k6NNlN  N(^2p^  —  k7NNz^NN2  +  k7NNlNN2  +  k8NN+NNz  —  k9N  N^2p^NN2  + 


kgNNNN2 


(3) 


d^(  5  —  ^v200^v2(a')  k3N n+NNz(A)  k5NN2^NN  k6NN2^NN  k7NNz^NN2 


(4) 


dN 


N2(a') 


dt  —  k4  NNz(A)N  N2^  k2  N  Nz  (a  ’)  N  (a  ■) 


d-Nfl 


—  k3NN+NN  —  k4NN+NN  —  k5NNl(A)NN  +  k5NNzNN  —  k6NN2^NN  —  k8NN+NN  +  k9NNlNN 


(5) 

(6) 


N(  ZP) 

~t  —  k6 NNzNN(  2 —  k9NN ^  2P)NN2 


(7) 


dN 


^=k9NN+NN2 


dN, 


dt 


■  —  —k3NN+NNz^  +  k4N  N+NNl  —  k8N  n+Nn 


(8) 

(9) 


dN 


dt 

dN 


3  —k3NN+NN  k4NN+NN 


NT 


kiNN+N„  +  k2NM+NR 


(10) 

(11) 
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dNe 

dt 


fciJV.r+AL  +  k7N  ,,+N„ 


(12) 


At  first,  simulation  of  Nitrogen  plasma  discharge  is  performed.  Here  the  chemical  kinetics  is  clearly  simplified 
and  Nitrogen  is  the  main  product  [11].  By  solving  the  steady  state  rate  balance  equations  we  calculate  the  species  density 
to  see  the  plasma  composition  on  atmospheric  pressure  and  at  temperature  300  K. 

To  get  the  solution  we  need  to  deduce  steady  state  balance  equations  from  the  differential  equations  in  previous 
section.  The  steady  state  balance  equations  are  derived  as 


^4  Nn+NNi  +  k5NNzNN  +  k6NNlN  N^2p^  —  k7NNz^NN2  +  k7NNzNNz  +  k8NN+NN2  k9N  N^2p^NNz 


+ 


k9NNNNz  =  0  (13) 

^n2(a)^ N2(a')  ~  1*3^ n+NNz(a)  —  —  k6NN2(A)NN  —  k7NNz^NNz  =  0  (14) 

— feliVW2(^)AfW2(a')  —  k2N  N2(a')N  w2(a)  =  ^  (1^) 

k3^N+NN  —  k4N  n+Nn  —  k5NNz^A')NN  +  k5NNzNN  —  —  k8N  n+Nn  +  k9NNzNN  =  0  (16) 

ks^N2  2p)  —  k9NN^  2P)NN2  —  0  (17) 


k,NN+NNz  =  0  (18) 

~k3N  n+NNz(A)  +  k4NN+NNz  —  k8N  n+Nn  =  0  (19) 

k3^N+NN  —  k4NN+NN  =  0  (20) 

K  NNtNe  +  k2NN+Ne  =  0  (21) 

k4NNtNe  +  k2NN+Ne  =  0  (22) 

3.  RESULTS 


The  Nitrogen  plasma  species  densities  from  pressure  1  torr  to  atmospheric  pressure  (760  torr)  at  300  K  are  shown 
in  figure  1.  All  the  species  density  for  various  pressures  was  determined  [11]  using  the  ideal  gas  law  i.  e.  PV  —  nRT , 
where,  P  is  the  pressure  of  the  gas,  V  is  the  volume  of  the  gas,n  is  the  amount  of  species  in  the  gas  mixture,  R  is  the  ideal, 
or  universal,  gas  constant,  equal  to  the  product  of  the  Boltzmann  constant  and  the  Avogadro  constant  and  T  is  the  absolute 
temperature  of  the  gas  [12].  From  the  result  we  see  that  at  1  torr  the  density  of  nitrogen  is  2.278  X  1016cm_3and  it 
increases  with  pressure  and  become  maximum  i.  e.  1.731  X  1019  cm~3  at  atmospheric  pressure.  Metastable  Nitrogen  and  atomic 
Nitrogen  densities  are  4.390  X  1015  cm" 3  and  1.723  X  1015  cm"3 at  1  torr  and  3.336  X  1018  cm"3  and  1.309  X  1018  cm"3 
at  atmospheric  pressure  respectively.  The  electron  density  is  5.504  X  1015cm_3at  1  torr  and  4.183  X  1018  cm_3at  atmospheric 
pressure.  From  the  figure  1,  it  is  observed  that  species  density  increase  rapidly  at  1  torr  to  50  torr.  After  50  torr  it  increases  slowly. 
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Figure  1:  Species  Density  Vs  Pressure  Curve  of  Nitrogen  plasma  at  300  K 

4.  CONCLUSIONS 

We  performed  a  mathematical  simulation  on  the  Nitrogen  plasma  discharge  to  find  the  species  for  various  pressures 
at  temperature  300  K.  The  species  density  of  various  Nitrogen  plasma  species  are  calculated  from  the  MATLAB  simulation. 
From  the  result  we  see  that  in  the  plasma  discharge  there  are  neutral  nitrogen,  atomic  nitrogen,  ionized  nitrogen  and  free 
electron.  The  electron  density  is  of  the  order  of  1015cm-3  at  1  torr  andl018cm~3at  atmospheric  pressure  and  at  300  K.  From 
the  figure  1  we  see  that  species  density  increase  rapidly  at  1  torr  to  50  torr.  After  50  ton-  it  increases  slowly.  It  is  observed  that 
a  significant  amount  of  Reactive  Nitrogen  Species  (RNS)  is  produced  in  the  Nitrogen  plasma  which  can  be  used  for  different 
applications. 
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